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Abstract 

We investigate the low-field phonon-limited mobility in armchair graphene nanoribbons (GNRs) 
using full-band electron and phonon dispersion relations. We show that lateral confinement sup- 
presses the intrinsic mobility of GNRs to values typical of common bulk semiconductors, and very 
far from the impressive experiments on 2D graphene. Suspended GNRs with a width of 1 nm 
exhibit a mobility close to 500 cm^/Vs at room temperature, whereas if the same GNRs are de- 
posited on Hf02 mobility is further reduced to about 60 cm^/Vs due to surface phonons. We also 
show the occurrence of polaron formation, leading to band gap renormalization of ~118 meV for 
1 nm-wide armchair GNRs. 
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Understanding the role of phonon scattering Qfl is of primary importance, since it 
provides information regarding the ultimate intrinsic mobility limit {fim) of a material, i.e., 
when all extrinsic scattering sources have been removed. This is especially important for new 
materials or nanostructured ones, such as graphene nanoribbons, where experiments are not 
fully comparable, given the presence of defects and other non idealities. Recent experiments 
have found /i^ in suspended graphene close to 10^ cm^/Vs near room temperature [5|, but 
with a sizable degradation after deposition on high-A; gate insulators 



coupling to the polar modes of the substrate pj-yj. 



3, probably due to 



Very few indications are available on mobility degradation in Graphene NanoRibbons 
(GNRs) deposited on different gate insulators |7|. As of now, /ij„ in sub-10 nm GNRs 
cannot be extracted from experiments since line-edge roughness (LER) is presently limiting 
mobility in state-of-the-art GNRs js- 10|. Theory, on the other hand, allows us to individually 
evaluate the impact of each scattering source on mobility, which is often a prohibitive task 
in experiments. Phonon- limited scattering is particularly important, since it represents the 
unavoidable scattering mechanism at finite temperature, and provides a good indication of 
the potential for the material in electronics, either for devices or interconnects. 

Here we investigate the effect of phonons and surface optical (SO) phonons on carrier 
transport in GNRs by means of a full band (FB) approach based on a tight-binding (TB) 



description 



11| of the electronic structure and of the phonon spectrum. We compute scatter- 



ing rates with first-order perturbation theory using the deformation potential approximation 
(DPA), and low- field mobility using the Kubo-Greenwood formula jl2|. We obtain the one- 
dimensional (ID) subbands of an armchair GNR from a TB p^-Hamiltonian accounting 
for energy relaxation at the edges [13:]. The electron energy dispersion is quantized in the 
transverse direction y with wavevectors fc^,^ = [2717]) /[{I + l)a], where a = 0.249 nm is the 
graphene lattice constant, / is the number of dimer lines and index t] runs from 1 to /. 
Graphene phonon spectrum is obtained with the force constant dynamic-matrix approach, 



including contributions up to the fourth nearest neighbors (4NNFC approach) ll| and us- 
ing force constant parameters extracted from first-principles calculations [l^. Each of the 
six phonon branches of graphene, labelled by the quantum number j, is splitted into / ID 
sub-branches, with transverse wave vector g^^ {qyp = (27r /?)/[(/ + l)a] for (3 = 0, ...,//2 — 1 
and qyp = [2Tt{/3 + 2)]/[{l + l)a] for /3 = 1/2, ...,/-!}. 

The momentum relaxation rate of an electron in the initial state k = [kx, kyrj) accounting 



for scattering from GNR phonons is obtained from the Fermi Golden Rule, summing over 



all final states k' 



(fc^, kyrj'), conserving total energy and longitudinal momentum 
1^ 1^1^ dq^ . ...J^R (1 + cos^kk' 
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where 



J — Dqp), if j is a longitudinal acoustic (in-plane optical) mode, 



Dac and Dqp = 1.4 x 10^^ eV/m |4| are the acoustic (AC) and optical (OP) deformation 
potentials, respectively, k'^ = kx ± Qx, = |q| = (qI + Qy/sV^'^, p ^ 7.6 x 10~^ g/cm^ is 
graphene mass per unit area [41, n~ is the Bose- Einstein occupation factor and riq = 
In addition, kp = 7r/(\/3a), W = a/2 (/ — 1) is the GNR width, (1 + cos 6'kk') is the spinor 
overlap 15|] and f{E) is the Fermi occupation factor. In Eq. ([T]) the upper sign is for phonon 
absorption (ABS) and the lower for phonon emission (EM). Gn,r]',i3 is the form factor due to 
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16|. 



the transverse momentum conservation uncertainty 

We have investigated the electrostatic coupling between electrons in the GNR channel 
and remote phonons of the substrate considering the GNR deposited on an oxide layer of 



thickness tox, width Wox and placed at a distance d = 0.4 nm [12|. The oxide is backgated 



le longitudinal 
21] for the two 



by an ideal metal. Since phonon modes are almost constant as a function of t 
SO phonon wavevector [3], we assume the same energies E^q as in Ref. 
considered SO phonon modes. In addition, electrons are confined in the plane, therefore the 
electron-SO phonon scattering rate reads [l2[: 



1 ' f+kp r 



LWoxe'^FjGe-'''^'' 



nl 6 (e^,-E^tEIo) il-k'x/kx' ^'^^^"^ 



l + cos6'kk' 



where Q = {Qx,Qy) represents the 2D SO phonon wavevector, Q = |Q|, the sum ruris 



over all SO phonon modes, F| oc 



-/ (LWox) is the electron-phonon coupling parameter |3l.ll2|. 
G{j],rj' ,Qy) is the form factor jl2| which reduces to Grf,n',i5 if Qy = %i3- 

In Eq. 02]), (Qx) is the GNR static dielectric function calculated within the Random 
Phase Approximation (RPA) in the size quantum limit j4, 12\. For 10 nm-wide GNRs, 
screening of the electric field due to polar vibrations is instead modeled by means of the 
2D RPA graphene static dielectric function 620 (Q) jl^ • From a numerical point of view, 
in order to reduce numerical noise, we have approximated 6 in Eqs. ([1]) and ([2]) with a 



Gaussian window of standard deviation AE and, for the lowest AC subbranches, a collisional 
broadening approach has been implemented considering = h/2 [1/t {k^qyp)]. Finally, 
we have computed the low-field mobihty by means of the Kubo-Greenwood formula 
accounting for ID transport [12I. 

DPA formally leads to a zero coupling with the transversal acoustic (TA) and fiexu- 
ral (ZA) phonon modes, so that only scattering with longitudinal acoustic (LA) modes is 
typically considered Q]. Theory 18 1 and Raman spectroscopy |l9^ have shown that ZA 



modes are negligible down to 130 K. However, classical results based on a TB description of 



electron-phonon coupling 201] and recent ab-initio calculations 2l| have demonstrated that 
TA modes play a comparable role as that of LA modes in degrading /ij„. On the other hand, 
a physical description of graphene taking into account long-range interaction between car- 
bon atoms highlights an off-diagonal coupling to the TA modes through the modulation of 
the hopping parameters, which is smaller than the on-diagonal deformation potential contri- 
bution |22| . We choose to adopt a physically consistent approach and use DPA considering 
electron coupling only with LA, LO and TO modes, rather than euristically reintroduce the 
contribution of TA modes. We use Dac = 10.9 eV, extracted from DFT calculation for the 
GNR family 3Z + 1 [23], rather than fitting experiments which are affected by uncontrolled 
mechanisms and actually lead to a large spread of the considered values for Dac R S 21 |. 



Low-field mobility is shown in Fig. [T^ as a function of the electron density for different 
widths, close to 500 cm^/Vs is found for 1 nm-wide GNR, exceeding by almost one 
order of magnitude the experimental mobility of GNRs j^, and the intrinsic phonon- 
limited mobility of silicon nanowires {24] of comparable size. We find that is mainly 
limited by backward scattering involving AC phonons, due to the large mode-dependent OP 
energy offset (~ 130-160 meV) [l^. Unlike in 2D graphene, where oc l/n2D 0], the 
lateral confinement in GNRs leads to a non-monotonic n2D-dependence as also observed in 
CNTs [25| (Fig. [1^). For small W , Hm increases with n2D, due to the reduction of final 
states available for scattering. For wider GNRs biased in the inversion regime, electrons can 
populate excited subbands opening additional channels for scattering, thus reducing 

In Fig. [2^, is plotted as a function of temperature T. Similarly to what has been 
observed in small-diameter CNTs |26|, in narrow GNRs the dependence on W and T can 
be expressed by means of the empirical relation jiiniW^T) = /io (300 K/T) (VT/l nm)°''^c' 
where /xq ~ 391 cm^/Vs and ct^c = 2.65, which is close to cx expected for narrow 



GNRs, since oc t/(W x DOS) oc G~^DOS~'^ and G oc l/W'^ and the density of states 
DOS oc l/y/W. Of course, for large W saturates to that of 2D graphene. Since AC 
in-plane phonons scattering is dominant and ~ kT/hu for kT ^ hu, is inversely 
proportional to T (Fig. The mean free path in the first subband (L^) is shown in Fig. [T}? 
as a function of and in Fig. I^b as a function of T, where (Lk) = {v (k) r (k)), (k) is 
the group velocity and the expectation value (•) has been computed in the Brillouin zone, 



considering /(I — /) as the distribution function [l6[. At T = 300 K, (Lk) is of the order of 
few fim for larger GNRs, as expected in graphene flakes, while it is ~ 10 nm for narrower 
GNRs (Figs. [2b). In addition, (Lk) oc 1/T, as (Figs. lib). 

The SO phono n-limited mobility function of n2D is shown in Figs. |3^-b for 

W smaller than 10 nm, considering GNRs deposited both on Si02 and on Hf02. As in 
graphene ^,^\, the higher the dielectric constant, the larger the mobility suppression due to 
SO phonon scattering. In particular, we observe fiex down to 700 cm^/Vs for Si02 (Fig. |3^) 
and 60 cm^/Vs for IIf02 (Fig. Eb), due to the smaller energy offset of the emission processes. 
As in CNTs [27], ^ex oc W^^^ with aso (~ 1.4-1.6) dependent on and smaller than a^c- 
For W smaller than 5 nm figx increases with electron concentration due to the impact of 
screening, whereas for W = 10.10 nm for higher concentrations n2D the increase of available 
modes for scattering reduces mobility. 

Comparing Figs. [S^-b with Fig. [T^, it can be observed that SO phonons play a sec- 
ondary role for very narrow GNRs on Si02 but they become predominant with increasing 
W > 2.5 nm roughly for n2D < 10^^ cm^^, whereas they are predominant for all n2D densities 
in GNRs on IIf02. Comparison with experiments shows that fi^x is larger by up to one order 
of magnitude than mobility measured on GNRs deposited on Si02 js] and by a factor three 
than mobility measured on 10-nm-wide GNRs integrated with ultrathin Hf02 dielectric 
This gives a rough estimation of the increase in mobility that could be achieved through 
fabrication technology improvements capable to suppress the present dominant scattering 
mechanisms (e.g. LER). As can be noted in Fig. ^ and as it also occurs in CNTs deposited 
on polar dielectrics [27], ^ex oc 1/T'^. In particular, for IIf02 7 ~ 3. Since oc 1/T, SO 
phonon scattering dominates transport roughly above 100 K for all W (Fig. [2b), as for CNTs 
on Si02 13. 

Finally, we focus on the polaronic energy shift SE]^^ due to the electron-phonon coupling, 



computed exploiting the second-order perturbation theory [l2|. Fig. ^ shows 6Ei^ as a 



function of i^k for the lowest two subbands for the W = 1.12 nm case. 5i?k is weakly 
energy dependent near the cutoff subband, is independent of T and increases sharply in 



Ts 



26| , mostly OP phonons contribute 



26|, the contribution to SEi^ from AC 



correspondence of intersubband transitions. As in CN^ 
to 6Ei^ (inset of Fig. H^). Instead, unlike in CNTs 
phonons exhibits few peaks due to the transverse momentum conservation uncertainty (inset 
of Fig. H^). The polaronic binding energy, i.e. the polaronic energy shift referred to the first 



conduction subband edge Eci, is 5Ei,=5E\^{E\^=Eci) and is almost 59 meV for 



GNRs, close to that obtained for semiconducting CNTs with the same number I [26| (Fig. Hh) 



nm-wide 



and corresponds to a band gap renormalization 2 SEi^ ^118 meV for 1-nm nanoribbons, and 
to a relative correction of -35% of the energy gap Eg of 10-nm nanoribbons (Fig. |3)d). 

In conclusion, we have proposed a very accurate full-band approach to evaluate low-field 
phonon-limited mobility in GNRs. We find that is close to 500 cm^/Vs in suspended 
1 nm-wide GNRs at room temperature, and is suppressed down to 60 cm^/Vs in 1-nm wide 
GNR deposited on Hf02, due to coupling with SO phonons. The result is important from 
the point of view of methodology and of fundamental physics, since the corresponding mean 
free paths range from 1 to 10 nm, undermining the possibility of performing ballistic or 
coherent transport experiments at non-cryogenic temperatures. It is also important from 
the application point of view: whereas suspended 2D graphene has an intrinsic mobility 
at room temperature several orders of magnitude larger than that of bulk semiconductors, 
narrow GNRs with reasonable semiconducting gap have only slightly larger mobility than 
comparable silicon nanowires. Finally, we also find polaron formation in armchair GNRs, 
with a remarkable band gap renormalization of up to 35% in the case of 10 nm-wide ribbons. 

This work was supported in part by the EC 7FP through NANOSIL (n. 216171), GRAND 
(n. 215752) grants, and by the MIUR-PRIN project GRANFET (Prot. 2008S2CLJ9). Au- 
thors thank www.nanohub.org for the provided computational resources. 
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FIG. 1. (a) and (b) (Lu) for an electron in the lowest subband as a function of for different 



W. DFT calculations 



2l[ | for graphene are also shown in (a). 
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FIG. 2. (a) mobility (b) mean free path (^k) for the lowest subband, and (c) mobility fi^x for 
GNR on Hf02 as a function of T for different W. In (a) and (c) dashed lines correspond to the 
empirical formula for fj,in{W,T) and solid lines in (c) to the fit for fj,ex- 
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FIG. 3. function of n2D for GNR deposited (a) on Si02 and (b) on Hf02. 
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FIG. 4. (a) Polaronic energy shift 5Ei^ as a function of the energy E]^ for the first two subbands 
(W = 1.12 nm). Sohd curves correspond to T= 300 K, dotted curves to T= 30 K. Inset: OP and 
AC (x 8) contributions to JE'k for the first subband. (b): Polaronic binding energy 6Eij (left) and 



polaronic correction to the band gap 25Eh/Eg (right) as a function of W. Results for CNTs 
are also reported. 
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